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Abstract 

In this paper we would like give a short review of our recent works and current 
investigation in modeling of solid film surfaces growth based on computer Monte Carlo 
simulations. We also discuss effect of some simplification in growth rules on speed- 
up of computations, and the time evolution of surface morphology (for example its 
roughness or anisotropy). We obtained qualitative agreement of computer simulation 
with experimental results for homoepitaxially grown simple metal films. 
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1 Introduction 

In the last 40 years history of surface growth science, computer simulations were par- 
ticularly useful and important for better understanding of modes and mechanisms of 
growth. Great number of simulations is based on so-called solid-on-solid (SOS) model 
approximation, when each particle seats on the top of the others and film surface may 
be described by a single- valued function of its height h{r) in planar substrate coordinate 
r (see [1-3] for review). In this paper, we present a short review of our recent works 
[4-6] as well as current investigation [7,8] and oncoming ones [9]. In all these papers we 
concentrate on films growth on lattice with simple cubic symmetry. 

Section 2 is devoted to presentation of the model. In Section 3 we present results 
which prove that a considerable speed-up of computations may be achieved at cost of 
some assumed simplification of growth rules and particle behaviors in computer sim- 
ulation. Then influence of model control parameters on surface morphology and its 
dynamics is presented in Section 3 and summarized in Section 4. 

2 Model 

Presented here results of simulations are based on simple SOS model of epitaxial thin 
films growth [4-8] based on Wolf- Villain model [10], in which after randomly chosen 
place of a particles initial contact with the surface, they relax to the site offering the 
strongest bonding (i.e. they maximize the number of particle-particle lateral bonds 
(PPLB)). The place of particle settling is probed among site of initial deposition plus 
their four nearest-neighbors. Then each particle repeats relaxation Ldif times and so 
Ldif may be considered as a range of diffusion. 
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This simple and deterministic rules of growth [10] are generalized in this paper by 
accounting on: 

• strengths of horizontal bonds {Jx, Jy) and vertical bonds {Sx, Sy) of PPLB in 
particle level layer and one layer below, respectively, 

• diffusion barrier {Vx, Vy), which restricts migration and stops it for infinitely large 
values of V, 

• temperature of substrate T. 

The temperature T (>iitries calculation by the Boltzmann factor exp{—E/kT) for 
probability p{E) of picking a site of energy E, which replaces the deterministic rule of 
fixing the particle at site of minimum E. Desorption process is neglected becouse typical 
desorption energies for metals are of the order of several eV. With such modification 
we may grow anisotropic surfaces for model control parameters {Sx ^ Sy, Jx Jy, 
and/or Vx 7^ Vy) [6,8]. The difi^usion barrier V reduces probability of movement by 
factor exp(14/fcT) or exp(Vj,/A:r) for considered direction x or y of migration. V values 
are always positive, while negative J and S correspond to the assumed tendency to 
maximize PPLB. Additionally, the presence of S bonds to particles below in particle 
energy may play role of the Ehrlich-Schwoebel barrier near the film edge [11,12]. The 
S term makes it difficult to reach the step edge and thus limiting horizontal growth of 
a cluster of atoms [5,13]. 

3 Results of Simulations 

The simulation were carried out on L x L large square lattice with periodic boundary 
conditions to minimize infiuence of boundary effects. The particles are represented as 
unit cubes which can occupy only discrete position in (2+l)-D space. 

3.1 Diffusion Range 

Firstly, we evaluate speed-up of computation when Monte Carlo simulations based 

on Arrhenius dynamics are reduced to oversimplified random deposition (RD) model 
through intermediate stages. If td is computing time to deposite a particle at the first 
contact site somewhere on the surface, then we need at least tuo — N ■ td time, where 
N is total number of the deposited particles, to complete simulation. The RD model 
correspond to Ldif = 0. Increasing Ldif = 1 results in changing model to local dif- 
fusion with single step relaxation and increase time of computation by the additional 
time tr necessary for relaxation procedure: tssR = N ■ td + N ■ tr- Case Ldif > 1 
corresponds to the local diffusion model when simulation time may be evaluated as 
tLD = N ■ td + N ■ Ldif ■ tr- Finally, for full-diffusion model {Ldif — > 00) we have 
tpD = N ■td + {1 + 2 + .. . + {N - 1) + N) - tr. Note, that almost deterministic rules of 
growth speed-up computation on a ratio k — tpo/tLD ~ {N + \)/2Ldif ~ OL^ /2Ldif , 
where the total number of particles A'^ — 8L^ is given by lattice lateral size L and film 
thickness/substrate coverage 6. And — as we shall show later — for such simplification 
we may loose some significant film properties and/or misinterpret particles dynamics. 

The range of diffusion Ldif influences surface roughness measured by standard sur- 
face width a, defined as the root-mean-square of surface heights. We found that 
increasing of diffusion range make surface smoother, however some saturation level 
Csat{Ldif — » 00) > is observed. The saturation level depends on strong of attrac- 
tive forces among particles which is governed by J [5]. The results are collected in 
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Table 1: Influence of bonding energy J on roughness saturation level a sat for large diffusion 
range Ldif 3> L. Lopt gives a rough estimate of the diffusion range for which saturation of a 
starts. = 10 ML, V = S = 0. 



J 


0.00 


-0.05 


-0.10 


-0.25 


-0.50 


-0.75 


-1.00 


-5.00 


-20.0 


Lopt 





1 


5 


«50 


wlOO 


W200 


w300 


w500 


> 5000 




3.16 


2.83 


2.53 


1.83 


1.30 


1.10 


1.00 


0.50 


0.36 



Tab. 1. Note, that effective difFusion range Le// may be also controlled by substrate 
temperature T as we show later in Section 3.4. 

3.2 Growth Dynamics 

The roughness a for film surfaces growing by deposition often obeys Family- Vicsek 
scaling law: a (x L°' ■ /(S/L^), whore scaling function is given by f(x) for x <C 1 

and f{x) 1 for a; ^ 1 [14,15]. For substrate low coverage ((? < 1 ML) also third regime 
is observed with growth exponent (3 = (3rd = 1/2, similar to simple RD [5,8,15,16]. 
Deviation from scaling hypotheses may be also observed when particles are allowed to 
climb on the top of other ones which results an unstable growth when surface width 
growth infinitely with increasing of coverage 6 [5,17,18]. This may bo considered as 
infinite value of roughness exponent (a oo). 

When also anisotropy in growth is included (e.g. for Ag/Ag(110) or Cu/Cu(110)), 
then we observe power law for roughness a oc 0^ for small 9 [8] (model parameters for 
Cu and Ag are taken after [19]). Before reaching 6 = 10* ML we still do not observe 
roughness saturation which may indicate that particles jumps to the higher levels may 
be not essential for presence of unstable growth {a —> oo) as it was earlier suggested 
in [5,18]. In both cases, either stable or unstable growth, surface roughness a increases 
initially as 9^. The growth exponent (3 may also depend on substrate crystallographic 
orientation as was shown by Elliot et al. for Ag(lll) and Ag(lOO) [20]. This can be 
explained only within models that allow for anisotropy. 

Family- Vicsek law fails for inhomogeneous multilayerd a/6/a-like sandwiches as a 
result of effective changes of growth rules — or more precisely — strength of interaction 
during deposition [7]. 

3.3 Surface Morphology Anisotropy 

As we mentioned in Section 2, the anisotropy may originate from either anisotropy in 
bonds {Sx 7^ Sy, Jx 7^ Jy) or anisotropy in barriers 14 7^ K/ [6,8]. Then anisotropy may 
be quantitatively measured from the difference in the height-height correlation function 
G(s) = {h{r + s) • h{r)) — {h{r))^ in two x and y directions, where (• • •) denotes spatial 
average over all sites r. Correlation G(1,0) are positive for terrace-like, smooth surfaces 
(in X- direction) while its negative values correspond to spiky and rough structures. Thus 
G{1,0) — G(0, 1) normalized to — G(0, 0) may be used as a quantitative measure of 
anisotropy e = [G(l, 0) — G(0, 1)]/G(0, 0) [6] and may be competitive to a direct pictures 
of the surface morphology obtained from simulations [8]. We found that anisotropy 
parameter e depends on the ratio of | Jx — J^l/I Ji; -|- Jj/| and not on the simple difference 
\Jx — Jy\ only. As it may be expected, the anisotropy vanishes for very large diffusion 
barriers Vx = Vy ^ 00, when RD results are reproduced. When the first monolayer 
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Table 2: Dependence of roughness a on substrate temperature T for Pt on Pt(lll) diffusion. 
Model control parameters are after Rcf. 13. 



T[K] 


300 


600 


900 


1200 


1500 


a [ML] 


2.04 


1.87 


1.74 


1.60 


1.49 



Table 3: Number M of mobile' particle's after L,/, / = ID"'' jumps for liomoepitaxial growth of 
Cu/Cu(110) and Ag/Ag(;il()). The mo(k-l i)ara.iiK-tcrs arc alter Rcf. 19. 



T [K] 


300 


600 


900 


1200 


1500 


Cu/Cu(110) 


1 


50 


239 


913 


1348 


Ag/Ag(110) 





28 


147 


529 


1069 



is completed, the anisotropy parameter e seems to become insensitive to increasing 
coverage [6]. 

3.4 Influence of Substrate Temperature on Film Properties 

We expect that increasing temperature should make the surface smoother since energetic 
barriers are easier to overcome and also surface bonds may be broken more easily. The 
effect of temperature increase is presented in Tab. 2. Model parameters are taken after 
[13] where simulations of Pt on Pt(lll) behaviors were studied. 

Dependence of effective diffusion range with substrate temperature may also bo 
easily predicted; for low temperatures particles freeze at the place of deposition and 
migration is stopped (T — > Lefj = 0). Then wc reproduce RD results with 
Poisson distribution of film heights. On contrary, hot substrate make particles mobile 
(T — » oo ^ Leff — » oo) and particles stop only if diffusion range is reached in the num- 
ber of subsequent jumps [8]. The quantity introduced to describe particle mobility is 
number M of particles which are still unsettled on the surface when the diffusion range 
is reached, that is the number M of particles which changed their position during the 
last step (Ldif) any iteration cycle. M is increasing function of substrate temperature 
T, as expected (see Tab. 3). Note, that the increase in Ldif plays similar role as an 
increase in substrate temperature T for stochastic deposition processes. 

Finally, we check how temperature influences surface morphology measured by e 
and applied to Cu/Cu(110) and Ag/Ag(110). For absolute zero temperature (T — » 0) 
we recover RD results (e — > 0). With increasing temperature the anisotropy coefficient 
reaches maximum and then drops, so that, we again have RD results e = for T ^ oo. 
In intermediate temperature surface morphology evolves from isolated small islands. 



Table 4: Influence of substrate temperature T on surface morphology anisotropy for 
Cu/Cu(110) and Ag/Ag(110) [8]. Model control parameters arc taken after Ref. 19. 



T[K] 


300 


600 


900 


1200 


1500 


s (Cu) 
£(Ag) 


0.275 
0.197 


0.732 
0.563 


0.828 
0.729 


0.535 
0.585 


0.508 
0.453 
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through strongly anisotropy long string of atoms along dimmer rows to larger and more 
isotropic islands also across dimmer rows [8,19]. 



4 Summary 

The described model proposed in [4-8] allows for qualitatively correct description of 
thin solid films grown by molecular beam epitaxy where dominant physical processes 
responsible for film growth is surface diffusion. We show that simplification of rules of 
growth — when next particle incoming only when position of a previous one was fixed — 
speeds up computation drastically particularly for large lattices and large coverage. This 
is justified for not very large flux of the incoming particles. Actually, bigger Ldif also 
mimics smaller flux of the deposited particles. The anisotropy in the model parameters 
allows to consider not only sc(lOO) symmetry surfaces, but also particles deposited on 
bcc(llO) surfaces [8]. 

From simulation we were able to reproduce data on thermal, time and structural 
evolution of surface morphology. Usually the surface roughness measured by standard 
deviation of film height increases with time (or coverage) accord to power law [5,8]. Such 
behavior is characteristic for early stages of growth of self-affine surfaces (see [21] for 
review) . 

We believe that it is important to study the fractal dimension for the (2-|-l)-D sur- 
faces for better understanding of the growth phenomena. We also intend to investigate 
more throughly the anisotropic growth [9]. 
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